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INTRODUCTION 


This  report  presents  an  evaluation  of  the  detection  performance 
for  the  time-varying  adaptation  minimum  power  adaptive  beamformer  (ABF). 
The  ABF  is  essentially  based  on  the  Lj  Norm  Least-Square  Formulation.  The 
ABF  algorithm  described  by  Shen  (1977)  was  modified  in  this  report  to  adapt 
to  various  levels  of  noise  in  the  operating  environment.  The  detection  per- 
formance study  is  based  on  the  modified  algorithm  in  which  the  design  param- 
eters are  kept  constant  to  demonstrate  the  algorithm's  ability  to  be  implemented 
for  operational  purposes. 

Another  element  discussed  in  this  report  is  the  problem  of  coda 
rejection.  The  ABF  capabilities  for  coda  rejection  can  be  automatically  ap- 
plied in  a front-end  detection  system  (Shen,  1978a,b).  The  coda  of  seismic 
signals  from  strong  earthquakes  consists  of  substantial  high  amplitude  and 
long  duration  wavetrains.  These  could  be  due  to  the  complex  nature  of  source 
functions  and/or  the  scatterings  of  waves  caused  by  the  complex  structure  in 
the  media  along  the  propagation  paths.  Under  these  circumstances,  an  auto- 
matic signal  detector  is  jammed,  or  continuously  yields  false  alarms  and 
loses  its  functions  for  a certain  period  of  time  before  returning  to  normal 
operation. 

Section  II  briefly  describes  the  new  ABF  algorithm.  Evaluation 
procedures  and  results  are  presented  in  Section  III.  The  data  base,  analysis, 
and  a statistical  estimate  of  the  detection  probability  curves  are  described  in 
detail  with  a summary  in  Subsection  III-C.  Conclusions  and  suggestions  are 
presented  in  Section  IV.  Section  V lists  the  references  cited  in  the  text  of  this 
report.  Finally,  Appendices  A and  B present  the  processed  data  displays. 
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SECTION  II 

THE  ADAPTIVE  BEAMFORMER 


A time-domain  linear  adaptive  beamformer  (ABF)  was  develop- 
ed to  be  used  in  an  operational  front-end  detection  system.  The  objective  is 
to  increase  the  detection  capability  of  a seismic  array  and  to  extract  the  low- 
magnitude  signals  of  unknown  waveform  from  noise.  The  algorithm  minimizes 
the  squared  filter  output  subject  to  the  constraints  which  pass  energy  from 
steered  location.  Mathematically,  let  y(t)  be  the  adaptive  filter  output  at 
time  t for  an  array  of  N sensors  (channels)  with  a 2J+1  points,  or  tap- 
pings of  filter  length  per  channel.  Then,  we  write 


N J 

y(t)  = £ 2 a(j)X(t-jAt) 

i=l  j=-J 


(II- 1) 


where  Xi  is  time-aligned  data  input  for  the  i channel.  At  is  the  sampling 

interval,  and  a.(j)  is  the  adaptive  filter  coefficient  for  the  i^1  channel  and 
th  ^ 

the  j point  of  the  channel  data.  In  the  algorithm,  a (j)  is  updated  each  data 
point  in  the  time -domain  by  the  expression 

a*+At(j)  = a*(j)  + £A(t)  y(t)  ^X(t-jAt)  - X.(t-jAt)J  (U-2) 


(II- 2) 


where  a.(j)  is  initialized  to  satisfy  the  Levin's  constraints  (Levin,  1964) 


fi  a'(i)  = V 

for  the  center  point  output  described  in  equation  (II- 1 ) and  where 


(II- 3) 


I 


1 1N 

X(t-jAt)  = — £ X.(t-jAt) 


(n-4) 


X is  introduced  to  preserve  the  constraints  which  were  initialized  at  the  very 
beginning  of  processing,  Ji  is  the  convergent  rate  and  is  an  input  parameter 
and,  finally  \(t)  is  a time-varying  and  adaptation  parameter.  The  algorithm 
was  formulated  (with  A(t)  = 1. 0)  by  Burg  et  al.  (1967)  and  has  been  investigated 
by  various  authors  (Kobayashi,  1970;  Frost,  1972;  Owsley,  1973;  and  Gangi  and 
Byun,  1976).  The  algorithm  achieved  substantial  noise  reduction  gain  relative 
to  the  beamsteer  (time-align  and  average)  processor.  However,  experimental 
work  on  the  recorded  seismic  data  indicated  a substantial  signal  loss  as  com- 
pared with  the  beamsteer  output.  That  resulted  in  the  optimum  signal-to- 
noise  ratio  (SNR)  gain  to  be  around  1 dB  relative  to  the  beamsteer  processor. 

A better  time-domain  adaptation  is  needed  to  achieve  a better 
SNR  gain.  It  was  found  that  a quasi  norm  adaptive  filter  improved  the  per- 
formance both  in  reducing  the  noise  and  in  preventing  the  desired  signal  from 
being  degraded.  Currently,  the  Lj  norm  ABF  has  been  modified  in  an  attempt 
to  make  the  ABF  processor  implementable  for  optimum  operation  in  various 
levels  of  seismic  noise.  Presently,  the  time-varying  adaptation  rate,  A(t) 
is  chosen  as  follows: 


x(t)  = yLTS(t) 

ySTS(t)  x p(t) 


(II- 5) 


where 


yLTS(t)  = (1_a)  yLTS(t_1)  +“ 


(H-6) 


ySTS(t)  = t1'/3)  ysTS(t_1)  ly(t)l 


(II- 7) 


a is  a long-term  smoothing  constant  of  filter  output  for  noise  adaptation  so 
that  /I  would  not  change  significantly  when  the  noise  level  changes  in  the  op- 
erating environment.  (3  is  the  short-term  smoothing  constant  of  filter  output 
intended  to  prevent  the  signal  loss.  It  is  noted  that  the  performance  of  the 
adaptive  beamformer  is  highly  dependent  on  a , (3,  and  /Z . 

Basically  an  norm  adaptive  filter  is  formulated  by  minimiz- 
ing  the  absolute  value  of  the  ABF  output.  This  is  in  contrast  with  the  more 
conventional  norm  adaptive  filter  formulated  by  minimizing  the  mean  square 
or  power  of  the  ABF  output.  According  to  Claerbout  (1976),  the  advantage  of 
an  Lj  norm  filter,  which  tends  to  find  the  median  of  data  samples,  is  that  such 
a solution  is  relatively  less  sensitive  to  large  deviations  of  the  data.  By 
weighting  the  time- varying  adaption  rate  A (t ) inversely  as  the  estimated 
signal-to-noise  ratio  of  the  ABF  output,  the  computation  of  the  update  of  all  of 
the  channel  filters  is  made  relatively  insensitive  to  large  fluctuations  in  the 
ABF  output.  Thus,  when  large  deviations  of  the  ABF  output  are  encountered 
due  to  signals  along  with  the  enhanced  signal  generated  noise,  the  filter  update 
is  automatically  suppressed  to  avoid  minimizing  the  power  of  the  signal  due  to 
the  enhanced  signal  generated  or  scattered  signal  component.  In  the  sense 
that  this  feedback  to  X(t)  makes  the  adaptive  filter  design  must  less  sensitive 
to  large  deviations  of  the  filter  output,  the  filter  design  is  considered  closer 
to  a minimum  than  a minimum  norm  optimum  solution.  Basically,  the 
formulation  is  a matched  filter  of  the  ABF  output  and  each  channel's  devia- 
tion from  the  conventional  beamform.  The  approximate  L.^  formulation  is  a 
matched  filter  between  the  sign  of  the  ABF  output  and  each  channel's  deviation 
from  a conventional  beamform.  Our  experience  indicates  that  this  largely 
eliminates  the  previous  problem  of  ABF  filters;  which  is  the  annihilation  of 
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SECTION  III 

DETECTION  PERFORMANCE 


l 


The  time-varying  adaptation  rate  adaptive  beamforming  algo- 
rithm was  evaluated  using  the  short-period  array  data  from  the  Korean  Seis- 
mic Research  Station  (KSRS).  The  KSRS  short- period  array  has  19  sensors 
with  an  aperture  of  about  10  km.  The  configuration  consists  of  two  concentric 
rings:  6 sensors  in  the  inner  ring  and  12  sensors  in  the  outer  ring  with  one  in 
the  center.  The  data  used  in  this  report  were  digitized  at  a sampling  interval 
At  = 0. 1 seconds. 


A.  CODA  REJECTION 


Figure  III-l  shows  the  beamforming  output  of  the  wide  band 
(0.  5-3.  5 Hz)  data.  The  signal  was  originated  from  the  Ceram  Sea  on  1 Jan- 
uary 1977.  The  first  trace  in  the  figure  is  the  beamsteer  signal  and  the  next 
two  traces  are  the  ABF  beams  with  convergent  rate  = 2.  0 and  jS  = 2.  5,  re- 
spectively. The  time-varying  adaptation  parameters  used  in  processing  are 
a = 0.  0045,  and  )3  = 0.  65  which  were  found,  by  various  tests,  to  have  opti- 
mum SNR  for  the  initial  P wave  arrival.  The  second  pulse  of  the  signal, 
which  arrived  11  seconds  after  the  initial  arrival,  is  the  pP  phase  of  the  P 
wave  signal  that  reflected  from  the  earth's  surface  near  the  source.  The 
secondary  phase  on  the  beamsteer  output  has  amplitudes  higher  than  the 
initial  phase,  suggesting  that  it  was  substantially  contaminated  with  the  multi- 
path  and  scattered  energy.  The  beams  of  the  next  two  traces  indicated  that 
the  contaminated  components  of  energy  were  amply  rejected  by  the  adaptive 
processor.  The  scattered  energy  after  the  pP  phase  was  also  substantially 
suppressed  in  the  ABF  beams  as  compared  with  the  beamsteer  output  in  the 
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figure.  This  is  due  to  the  fact  that  the  initial  arrival  of  the  seismic  wavetrains 
may  have  better  coherency  among  the  channels,  while  the  multi- scattered  en- 
ergy in  general  shows  poorer  coherency  and  therefore  was  rejected  by  the 
adaptive  beamformer. 

Figure  III- 2 shows  the  beamforming  outputs  for  the  same  signal 
shown  on  Figure  III- 1 with  input  data  in  the  0.  5-1. 5 Hz  passband.  Rejection 
of  the  coda  by  the  adaptive  filter  is  more  pronounced  in  the  narrower  band  data. 
The  pulse  which  arrived  1 minute  35  seconds  later  after  the  initial  P wave  ar- 
rival is  PP,  the  reflected  phase  of  the  signal  on  the  earth's  surface  approxi- 
mately at  the  mid-point  between  the  source  location  and  the  receiving  array. 
The  third  identifiable  pulse  on  the  ABF  beam  which  arrived  2 minutes  and  15 
seconds  later  after  the  primary  signal  is  possibly  the  PcP  phase  which  is  due 
to  the  reflection  from  the  core  of  the  earth. 

To  investigate  how  the  beamsteer  and  ABF  processors  respond 
to  the  wavetrains  of  different  coherency  among  the  channels,  the  data  from  the 
Ceram  earthquake  were  used  to  form  the  array  beam  in  various  look  azimuths 
with  a fixed  apparent  P wave  velocity.  Figure  III- 3 illustrates  the  beam  pat- 
tern of  the  measured  peak-to-peak  amplitude  for  the  first  signal  pulse  shown 

in  Figure  M-l  (in  the  0.  5-3.  5 Hz  passband).  The  decibels  are  computed  rel- 

o 

ative  to  the  beamsteer  output  at  the  182  azimuth.  Because  the  ABF  is  design- 
ed with  the  time- varying  adaptation  rate,  the  adaptive  filter  (fi  = 2.  0)  response 
in  general  follows  that  of  the  beamsteer  processor  in  this  figure.  However, 
the  measured  beam  patterns  on  the  basis  of  the  second  pulse  signal  in  Figure 
III-l  does  show  significant  difference  between  the  beamsteer  and  ABF  pro- 
cessors as  shown  in  Figure  HI-4.  That  is  why  the  adaptive  filter  is  able  to 
achieve  more  coda  rejection  than  the  beamsteer  filter.  The  response  pattern 
actually  measured  on  the  basis  of  recorded  data  could  slightly  differ  from  one 
signal  to  the  other.  The  earthquake  event  used  in  the  study  was  randomly 
selected  to  demonstrate  coda  rejection. 
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BEAMFORMING  OUTPUTS  FOR  A STRONG  EVENT  FROM  THE  CERAM  SEA 
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B.  DETECTION  PROBABILITY  MEASUREMENT 

The  main  task  in  this  report  is  to  measure  the  detection  prob- 
ability as  a function  of  bodywave  magnitude  (m^)  for  the  ABF.  This  involves 
a statistical  estimate  and  usually  requires  a substantial  data  base  to  form 
reliable  and  meaningful  results. 

1.  Data  Base 

A total  of  129  events  were  processed  and  used  in  this  evaluation. 
First,  events  listed  in  the  January  1977  monthly  listing  of  Preliminary  Deter- 
mination of  Epicenters  by  the  U.S.  G.  S.  National  Earthquake  Information 
Service  (NEIS)  were  selected,  with  the  following  criteria  applied: 

20°  < A < 100° 

where  A is  the  epicentral  distance  of  a source  location  to  KSRS.  This  constraint 

was  imposed  on  event  selection  because  of  the  intention  to  design  the  detector 

for  teleseismic  P waves.  The  processing  was  performed  on  95  events  from 

January  1 to  15,  1977.  Table  III  - 1 lists  the  event  information  from  the  NEIS 

bulletin  and  our  own  measurements  of  m,  . There  are  four  columns  on  Table 

b 

III-  1 concerning  the  m^.  From  the  left,  the  first  is  from  the  NEIS  bulletin, 
the  second  is  the  single-site  (SS)  measurement,  and  the  third  and  fourth  columns 
are  the  measurements  of  m by  the  conventional  beamsteer  (BS)  and  the  ABF, 

D 

respectively.  A blank  in  a row  means  no  detection.  However,  if  a SS  magnitude 
is  listed  (i.  e.  , detected),  array  beamforming  was  not  performed.  Therefore, 
neither  BS  nor  ABF  units  were  available  for  that  event.  But,  in  the  case  of 
detection  probability  estimates,  detections  were  claimed  for  both  the  BS  and 
the  ABF.  ID,  the  first  column  in  Table  III  - 1 , is  the  identification  number  which 
was  used  for  convenience  in  referring  to  the  processed  traces  in  Appendix  A 
and  Appendix  B. 
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PROCESSED  EVENT  LIST  FROM  NEIS  BULLETIN 
(01/01/77  to  01/15/77) 

(PAGE  1 OF  4) 
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PROCESSED  EVENT  LIST  FROM  NEIS  BULLETIN 
(01/01/77  to  01/15/77) 

(PAGE  3 OF  4) 
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Next,  a number  of  events  were  selected  from  the  January  1978 
NORSAR  bulletin.  The  purpose  was  to  select  some  low-magnitude  events  for 
the  ABF  evaluation.  Three  constraints  were  imposed  on  the  event  selection 
as  follows: 

<1>  NORSAR  mb  < 4.5 

<2>  NORSAR  Bulletin  V < 14.0  km/sec 

P 

<3>  Event-KSRS  Distance  20°  < A < 99° 

where  is  the  apparent  P wave  velocity  at  NORSAR  and  was  listed  in  the 
NORSAR  event  bulletin.  The  first  constraint  was  intended  to  improve  the 
ABF  evaluation  by  selecting  additional  low-magnitude  events  not  available  on 
the  NEIS  list.  The  second  constraint  was  applied  to  ensure  a satisfactory  de- 
gree of  reliability  for  a given  event  listed  in  the  bulletin.  Finally,  the  third 
constraint  was  imposed  to  limit  the  evaluation  to  teleseismic  P waves.  Pro- 
cessing was  performed  on  34  events  from  January  1 to  31,  1978.  Table  III-2 
lists  the  information  from  the  NORSAR  bulletin  and  the  measured  m . Appen- 
dix  B presents  the  processed  traces  for  those  events.  (Note:  Event  ID  was 
introduced  in  Table  HI-2  and  Appendix  B for  the  convenience  of  the  reader.) 

2.  Comparison  of  Station  m and  Bulletin  m, 

b b 

Measurement  of  bodywave  m^  magnitude  was  done  when  a signal 
was  detected.  An  analyst-type  of  detection  criteria  was  used  in  which  signal 
amplitude  rise  was  about  twice  that  of  the  maximum  noise  amplitude,  in  a one 
minute  time  gate.  The  measured  m^  was  computed  by  the  expression 

mfe  = log  + p (A,  h) 

where  A is  the  peak-to-peak  amplitude  in  m fi,  T is  the  period  and 

P (A,  h)  is  a correction  factor,  which  is  a function  of  epicentral 
distance,  A (in  degree)  and  the  depth,  h (in  km). 

The  P factors  in  Table  2 of  Veith  and  Clawson  (1972)  were  used  for  the  com- 
putations. 
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The  KSRS  short-period  array  has  an  instrument  response  peak 
centered  at  1.0Hz  (0.488  m/x /computer  count)  and  thereby,  T = 1.0  second 
was  used  for  all  nr^  computations.  The  NEIS  bulletin  lists  depths  while  the 
NORSAR  bulletin  does  not.  In  the  absence  of  depth  from  the  NORSAR  bulletin, 
the  40  km  'normal  depth'  was  used  for  P (A,  h)  corrections. 

Figure  III- 5 compares  the  site  1 measured  m,  at  the  KSRS 

D 

with  the  NEIS  or  NORSAR  m (numerical  m,  units  are  listed  in  Tables  III -1 

b b 

and  III-2,  and  are  separately  and  repeatedly  listed  in  the  seismic  trace  figures 

in  Appendices  A and  B for  convenience  to  the  reader).  The  dashed  slant  line  in 

Figure  III-5  indicates  the  difference  between  the  average  bulletin  m^  and  the 

average  measured  m,  at  site  1.  In  other  words,  the  measured  m units  at 

D b 

site  1,  on  the  average,  yielded  0.  35  units  higher  than  those  in  the  bulletins. 

A number  of  factors  could  possibly  attribute  to  the  bias  or  difference  in  the 
averaged  units.  First,  single  seismometer  measurement  tends  to  in- 

crease the  m^  unit  by  noise  contribution  to  the  signal  amplitude  (Note:  meas- 
urement was  done  by  peak-to-peak  amplitudes).  Second,  NORSAR  body- 
wave  magnitude  in  general  was  about  0.  2-0.  3 unit  smaller  than  the  NEIS  m 

b 

(Ringdal  and  Whitelaw,  1973). 

Figure  in  - 6 shows  the  comparison  of  NEIS  or  NORSAR  m 

b 

versus  the  measured  m by  the  conventional  beamsteer  at  the  KSRS.  The 

b 

slant  line  is 

Beamsteer  mfc  = -0.1+  NEIS  or  NORSAR  m^ 

which  indicated  that  the  average  KSRS  beamsteer  m^  was  0.  1 unit  smaller 
than  those  of  the  bulletin.  This  difference  is  due  to  the  beamforming  loss  at 
the  KSRS  short-period  array.  The  same  method  of  comparing  the  bulletin 
with  ABF  m^  is  shown  in  Figure  III -7  where  the  dashed  slant  line  is 

ABF  mb  = - 0.  1 + NEIS  or  NORSAR  ny 

Considering  the  average  of  all  signals  detected  by  the  ABF,  the  ABF  was 
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0.  1 unit  smaller  than  those  in  the  bulletin.  This  difference  is  the  same  mag- 
nitude of  loss  as  the  conventional  beamsteer. 

3.  ABF,  Beamsteer,  and  Single  Site  Detection  Probabilities  at 
KSRS 

Histograms  of  detection  and  non-detection  as  a function  of  body- 
wave  magnitudes  m^  were  generated  from  the  results  listed  in  Tables  III-l 
and  III- 2.  For  a given  event  in  the  tables,  if  a signal  was  detected  by  site 
l(i.  e.  , m^  value  is  printed  in  SS  column),  detections  are  automatically 

claimed  for  the  beamsteer  and  ABF  processors.  The  detection  percentages 
(in  terms  of  m^)  were  fitted  to  Gaussian  probabilities  by  use  of  a least- 
square  procedure  with  uniform  unity  correct  decision  probability  (Shen,  1976). 
Curves  of  90%  confidence  were  also  computed. 

Figure  III-8  shows  the  KSRS  site  1 detection  and  non-detection 
histogram  and  detection  probability.  In  the  figure,  MB50  means,  in  terms  of 
mb’  t*lC  ^etect^on  probability,  and  MB90  the  90%  detection  probability. 

In  this  case,  the  bodywave  magnitude  for  a 50%  detection  probability  is  a mean 
value  fimb  = 5.23,  with  a standard  deviation  crm^=  0.34. 

For  the  array  beamforming  results,  Figure  III-9  shows  the 

KSRS  beamsteer  detection  and  non-detection  histogram  and  detection  probability. 

The  50%  detection  probability  is  /Um  = 4.  62  and  the  standard  deviation  is 

b 

(7  = 0.  58.  The  conventional  beamforming  improvement  at  the  KSRS  short- 

period  array  is  0.  61  m^  units  or  12.  2 dB  for  the  MB50  and  0.  30  units 

or  6 dB  for  the  MB90.  The  12.  2 dB  beamforming  detection  improvement 
agrees  with  >/n"  array  gain,  where  N is  the  number  of  channels  of  an  array. 

(i.  e.  , 12.78  dB  for  N = 19). 

The  array  detection  capability  was  enhanced  further  by  the  ABF. 
Figure  III- 10  shows  KSRS  ABF  detection  and  non- detection  histogram  and  de- 
tection probability  where  MB50  is  4.  14  and  MB90  4.  84,  with  orm^  0*  54. 
Detection  improvement  of  the  ABF  over  the  conventional  beamsteer  is  0.  48. 
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unite  or  9.  6 dB  at  MB50  and  0.  53  units,  or  10.  6 dB  at  MB90.  This 
demonstrated  ABF  improvements  enhanced  the  KSRS  (19  sites)  detection  ca- 
pability to  the  equivalence  of  a conventional  beamsteer  operating  on  an  array 
consisting  of  190  sensors.  That  increase  in  capability  is  significant. 


It  should  be  noted  that  the  detection  improvement  (in  dB  or  m 

b 

units)  mentioned  here  was  the  threshold  reduction  rather  than  SNR  processing 
gain.  One  would  expect  that  the  SNR  processing  gain  should  be  equal  to  the 
threshold  reduction  if  processing  is  linear,  in  terms  of  input-output  SNR  re- 
lationship, such  as  the  conventional  beamsteer.  Previous  study  by  adding 
scaled  noise  to  signal  (Shen,  1977)  indicated  that  the  ABF  processing  gain  was 
not  equal  to  its  threshold  reduction. 


Bodywave  magnitudes  used  in  Figures  III- 8,  III- 9,  and  El- 10 
were  directly  taken  from  NEIS  or  NORSAR  bulletins.  The  events  vary  over 
a wide  range  of  epicentral  distances  and  thus  potentially  represent  a wide 
range  of  signal  amplitudes  corresponding  to  each  m^  value.  Because  the  epi- 
central distances  range  from  20°  to  99°,  it  was  felt  that  the  detection  proba- 
bility might  not  be  well-behaved  in  terms  of  Gaussian  distribution  and  result 
in  a large  event  population  standard  deviation.  Therefore,  for  the  purpose  of 
statistical  study,  bodywave  magnitude  was  adjusted  to  an  equivalent  50°  epi- 
central distance  from  the  KSRS  by  the  transformation 

mb(50,  h)  = mb(A,  h)  - p(A,  h)  + p(50,  h) 

where  A is  epicentral  delta  (in  degrees)  and  h is  depth  (in  km),  m (A,  h)  is 

b 

the  bulletin  magnitude,  and  p(A,  h)  is  a correction  factor  for  A and  h . 

Figure  III- 11  shows  the  site  1 detection  and  non- detection  his- 
togram and  detection  probability  where  the  magnitudes,  have  been  adjust- 
ed for  the  location  being  at  50°  from  the  KSRS.  The  adjusted  Gaussian  param- 
= 5.  18  and  crmb  = 0.  29.  The  same  procedure  was  accomplish- 
ed for  the  beamsteer  and  the  ABF.  Figure  III- 12  shows  the  histogram  and 
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detection  probability  for  the  beamsteer  results.  The  standard  deviation  was 

reduced  somewhat  from  that  of  Figure  III- 9,  but  only  by  0.  07  units.  For  the 

ABF  detection  results,  the  same  illustration  is  shown  in  Figure  III-13  where 

the  least-square  fitted  Gaussian  parameters  are  fim  = 4.05,  and  o"m  = 0.47. 

b b 

The  apparently  large  standard  deviations  of  the  conventional  beamsteer  and 
ABF  event  populations  are  probably  due  to  the  manner  in  which  NEIS  and 
NORSAR  m^  measurements  were  mixed.  Apparent  negative  magnitude  anom- 
alies of  NORSAR  events  of  4.  5 would  result  in  a stretching  of  the  prob- 

ability of  detection  curves  for  events  smaller  than  4.  5.  This  effect  did  not 
occur  on  the  single- site  curve  because  only  one  small  event  was  detected. 

C.  SUMMARY 

The  detection  performance  of  the  L,j  norm  adaptive  beamform- 
er  has  been  evaluated  using  the  KSRS  short-period  array  data.  It  has  been 
demonstrated  that  the  ABF  is  able  to  suppress  the  coda,  the  scattered  energy 
after  the  initial  arrival  of  a P wave,  much  better  than  the  conventional  beam- 
steer. In  the  0.  5-1. 5 Hz  passband,  the  signal-to-noise  ratio  gain  relative  to 
the  beamsteer  in  general  is  above  10  dB.  Also,  one  of  the  important  achieve- 
ments of  the  ABF  is  that  it  is  able  to  suppress  scattered  wavetrains  while  en- 
hancing the  secondary  P wave  phases  such  as  PP,  PcP,  and  pP. 

Bodywave  magnitude  m,  was  also  measured  on  the  basis  of 

b 

peak-to-peak  amplitude  with  correction  for  source  epicentral  distance  and 

depth.  On  the  average,  the  single  site  KSRS  m,  yielded  0.  35  m,  units  great- 

b b 

er  than  the  NEIS  and  NORSAR  m^.  The  beamsteer  processor  and  the  ABF 
both  yielded  0.  1 units  less  than  NEIS  and  NORSAR. 

On  the  detection  probability  measurement  evaluation,  a total 
of  129  events  recorded  at  the  KSRS  were  used  for  a detection  performance 
study.  Measurements  of  the  detection  capability  were  accomplished  for  a 
single- site,  the  conventional  beamsteer,  and  the  ABF.  Table  III- 3 summarizes 
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the  detection  probability  estimates  by  fitting  the  detection  percentages  of 
to  the  Gaussian  probability.  The  table  lists  50%  detectability  (mean  value, 
fi)  90%  detectability  m^,  and  the  event  population  standard  deviation 
SNR  improvement  of  the  conventional  array  beam  over  the  single  site  and  the 
ABF  over  the  conventional  beamsteer  were  about  12  dB  and  10  dB,  respective- 
ly. The  10  dB  ABF  detection  improvement  is  equivalent  to  a conventional 
beamsteer  SNR  gain  of  an  array  with  about  10  times  the  number  of  sensors. 

Finally,  the  data  used  in  this  report  were  recorded  in  the  winter 
season.  The  KSRS  array  may  have  better  performance  in  the  summer  season 
when  the  noise  is  quieter. 

These  results  are  based  on  subjective  analyst  interpretation  of 
the  detection  status  of  beamformed  and  adaptive  beamformed  processed  events. 
More  data,  with  automatic  detection  criteria  used  to  determine  the  detection 
status  of  events,  are  needed  to  more  accurately  estimate  the  detection  gains 
possible  by  applying  the  ABF  as  a front-end  seismic  event  detector. 


TABLE  III- 3 

SUMMARY  OF  DETECTION  PROBABILITIES  AND  IMPROVEMENTS 
FOR  ABF,  BEAMSTEER,  AND  SINGLE  SITE 


50% 

Detectability 
('im  b) 

90% 

Detectability 
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Deviation 

(aN 

Comments 

Site  1 

5.  23 

5.  67 

0.  34 

Bulletin 

Beamsteer 

4.  62 

5.  37 

0.  58 

"'b 

ABF 

4.  14 

4.  84 

0.  54 

Beamsteer 

Improvement 

Over  Single  Site 

0.61 

(12.  2 dB) 

0.  30 

(0.60  dB) 
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(9.  6 dB) 

0.  53 
(6.0  dB) 
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5.  18 

5.  56 

0.  29 

Bulletin 

Beamsteer 

4.53 

5.  19 

0.  51 

m^  are 

ABF 

4.05 

4.  65 

0.47 
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Beamsteer 

Improvement 

Over  Single  Site 

0.65 

(13.0  dB) 
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(7.4  dB) 

— 
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SECTION  IV 

CONCLUSION  AND  SUGGESTIONS 


A.  CONCLUSION 

The  ABF  was  developed  as  an  array  processor  for  signal  de- 
tection and  extraction.  There  is  no  need  to  inject  a pilot  signal  or  noise  sta- 
tistics to  train  the  system.  The  ABF  can  be  implemented  as  a front-end  de- 
tector. 

Using  the  KSRS  short-period  array  data,  evaluation  for  its  de- 
tection performance  has  been  accomplished  in  terms  of  its  capability  to  sup- 
press coda  wavetrains  or  the  scattered  energy  following  the  initial  arrival, 
to  enhance  the  secondary  P wave  phases,  accuracy  of  its  measurement, 
and  estimates  of  its  detection  probability  as  compared  with  the  conventional 
beamsteer.  A total  of  129  events  were  used.  On  the  average,  the  ABF  yield- 
ed the  same  magnitude  difference  of  0.  1 m units  as  the  beamsteer  when  com- 

b 

pared  with  the  bulletin  m^.  However,  the  ABF  improved  the  KSRS  detection 
capability  by  0.  5 m^  units  from  the  conventional  beamsteer  as  suggested  from 
the  statistical  estimate  of  the  processed  results. 

B.  SUGGESTIONS  FOR  FURTHER  STUDY 

The  present  ABF  algorithm,  programmed  for  and  catalogued 
in  the  IBM/360  TS44  system,  was  used  in  the  real-time  mode  except  for  off- 
line processing  using  recorded  field  tapes  and  was  limited  to  one  beam  (i.  e.  , 
one  azimuth  and  one  velocity).  The  processed  data  were  stored  in  a data  buf- 
fer for  a Calcomp  plot  and  for  other  purposes.  The  ABF  can  be  improved  to 
speed-up  computation.  To  implement  the  ABF  in  a real-time  front-end 

i 
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detection  system,  additional  software  work  is  needed.  Examples  of  this  soft- 
ware work  are  an  automatic  control  of  data  quality  among  the  channels  in  an 
array,  power  detector  decision  rules,  and  multiple  beams  to  cover  a certain 


area. 


The  Lj  norm  ABF  has  proved  very  successful  for  a seismic 
short-period  array.  On  the  basis  of  short-period  and  long-period  array 
studies,  it  is  expected  that  the  present  ABF  would  have  better  performance 
for  the  long- period  array  than  for  the  short-period  array.  Application  of  the 
Lj  norm  ABF  to  long-period  seismic  arrays  for  surface  wave  detection  and 
extraction  would  enhance  the  detection  and  discrimination  capability  of  seis- 
mic arrays. 

The  recent  growth  of  interest  in  regional  seismic  signal  detec- 
tion and  identification  strongly  suggests  that  applications  of  the  ABF  in  small 
arrays  (5  to  7 channels)  for  both  long-period  and  short-period  data  may  pro- 
vide an  effective  and  economical  solution  to  the  problem.  A feasibility  study 
and  investigation  of  this  application  would  be  beneficial. 
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PROCESSED  DATA  PLOTS  FROM  NEIS  BULLETIN 
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This  appendix  presents  the  processed  data  plots  for  those 
events  taken  from  the  NEIS  bulletin  dated  January  1 to  15,  1977.  Time  in  the 
figure  is  depicted  as  Julian  date,  hour,  minutes,  and  seconds.  In  general, 
the  computed  P-wave  arrival  is  70  seconds  after  the  start  time.  In  the  pre- 
sent size  of  figure,  it  is  about  10  seconds  for  0.  5 inch  scale.  The  line  of 
computer  print  on  the  top  of  each  figure  also  appears  in  Table  HI-1  herein. 
The  parameters  a = 0.  003  and  /3  = 0.  5 were  used  for  processing. 
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APPENDIX  B 

PROCESSED  DATA  PLOTS  FROM  NORSAR  BULLETIN 

| The  second  appendix  shows  the  processed  data  plots  for  the 

events  taken  from  NORSAR  bulletin  dated  January  1 to  31,  1978.  The  selec- 
| tion  criteria  and  brief  explanation  of  the  plots  have  been  discussed  in  an  earl- 

ier portion  of  this  report. 
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